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\^ \ A first measurement of tiie di-liadron interference fragmentation function of light 

quarks in pion pairs with the Belle detector is presented. The chiral odd nature 
of this fragmentation function allows the use as a quark polarimeter sensitive to 
1) \ the transverse polarization of the fragmenting quark. Therefore it can be used 

Q^' together with data taken at fixed target and collider experiments to extract the 

^ . quark transversity distribution. A sample consisting of 711 x 10^ di-hadron pairs 

was extracted from 661 fh^^ of data recorded near the T(4S) resonance delivered 
by the KEKB e+e" colhder. 
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^ \ The di-hadron interference fragmentation function (IFF), suggested first by Collins, Hep- 

pelmann and Ladinsky [1] describes the production of unpolarized hadron pairs in a jet 
from a transversely polarized quark. The transverse polarization is translated into an 
azimuthal modulation of the yields of hadron pairs around the jet axis. In addition the 
5^ \ IFF is chiral odd and can therefore act as a partner for the likewise chiral odd quark 

transversity function. The resulting amplitude is chiral even and therefore leads to ob- 
servable effects in semi deep inclusive scattering off a transversely polarized target [2,3] 
or in proton-proton collisions [4] in which one beam is transversely polarized. Besides the 
fragmentation into two hadrons, the fragmentation of a transversely polarized quark into 
one unpolarized hadron can be used to extract transversity via the Collins effect. Mea- 
surement of this effect at Belle [5] made the first extraction of transversity possible [6]. 
However as compared with the Collins fragmentation function, using the IFF to extract 
transversity exhibits a number of advantages. These are connected to the additional de- 
gree of freedom provided by the second hadron. It allows to define the azimuthal angle 
between the two hadrons as an observable in the transverse plane and at the same time 
integrate over transverse momenta of the quarks and hadrons involved. Because trans- 
verse momenta are integrated over collinear schemes in factorization and evolution can 
be used which are known and which do not need assumptions of the intrinsic transverse 
momenta [7]. Since the IFF is not a transverse momentum dependent function (TMD) 
it is universal and therefore directly applicable to SIDIS and proton proton data. From 
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both types of experiment results are available [2-4]. The results from SIDIS are indicating 
a non-zero IFF while the first analysis of the IFF effect at PHENIX [4] opens a way to 
disentangle transverse spin effects in proton proton collisions. Because the IFF is not a 
TMD function, extraction of transversity times IFF is not dependent on a model of the 
transverse momentum dependence, which is the case in the measurement of the Collins 
effect [6]. Here transverse momentum in the final state originates from a convolution of 
quark distribution and fragmentation function. This leads also to the Sudakov suppres- 
sion of the effect [8]. Technically, the extraction of the IFF from electron colliders is 
also easier, as the signal is not competing with asymmetries from QCD radiation. Also 
acceptance effects are smaller for the relative azimuthal angles of hadron pairs. 

2 Observables in e^e~ collisions 

The transverse polarization of the fragmenting quark leads to a cosine modulation of the 
azimuthal angle of the plane spanned by the two hadrons hi, h2 which is described by the 
vector R = P^ji — lying in the plane. Since the electron beams are unpolarized any 
effect that one would measure in one hemisphere of an event would average out. Instead 
one can make use of the fact that the spins of quark and anti quark in electron-positron 
annihilation are 100% correlated. Thus the correlation of the azimuthal angles of the 
vectors R° in the hemispheres a G {1,2} around the thrust axis with regard to the event 
plane is sensitive to the IFF. The measurement uses the center of mass system and defines 
the event plane as the plane which contains the beam axis z and thrust axis n. Figure 
[1] shows the coordinate system with the relevant quantities. With these the angles (pa 
between and the event plane can be expressed as 




Figure 1. Azimuthal angle definition. Azimuthal angles (pi and (p2 are defined relative to the 
thrust axis. 



0{i,2} = sgn [n ■ (z X n X (n X Ri,2)}] 

/ z X n n X Ri.2 \ 

X arccos rr ■ ^ ' . . (1) 

\|z X n| |n X Ri,2|y 

The product of the quark and anti quark interference fragmentation functions Hi ■ Hi is 
then proportional to the amplitude ai2 of the modulation cos(0i -t- ^2) of the di-hadron 
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pair yields [9]. The di-hadron fragmentation functions are dependent on the kinematic 
variables minv, the invariant mass of the hadron pair, and z = the normalized energy 
of the hadron pair. Here Eh is the energy of the hadron and Q the absolute energy 
transferred by the virtual photon. Therefore the normalized azimuthal yield of di-hadron 
pairs can be described as 

N{y, zi, Z2, mi, m2, 0i + ^2) oc a^iy, zi, Z2, mi, 7712) ■ cos(0i + ^2) (2) 

with 

auiy, zi, Z2, mi, m2) ~ B{y) ■ -. (3) 

e^gDi{zi,mi)Di{z2, mg) 

Here the sum goes over all quark flavors q and B{y) = i+cos^ e kinematic 

2 y^y 

factor describing the transverse polarization of the quark-anti quark with regard to its 
momentum. The polar angle 6 is defined between the electron axis and the thrust axis as 
shown in fig. [H The yields are dependent on zi, Z2, mi, m2, the fractional energies and 
invariant masses of the hadron pairs in the first and second hemisphere, and y. In the 
following the yields are integrated over y in the limits of the fiducial cuts on 6 introduced 
in section 12. 3[ The labeling of the hemisphere is at random, and it is experimentally not 
possible to distinguish between quark and anti quark fragmentation. Further information 
about the IFF can be learned from the dependence on the decay angle dh in the CMS of 
the two hadrons produced. Using a partial wave decomposition to isolate components that 
contain the interference between waves with one unit difference in angular momentum, one 
expects a dependence on sin 6h which gives the interference term between s and p waves. 
This term is expected to dominate at Belle kinematics and is favored by the acceptance. 
Since the acceptance is symmetric around 9h = j the p-p contribution proportional to 
cos6h should average out. Table [3] shows that this is approximately the case. 



2.1 Models 

As described in the previous section the IFF is an interference effect between hadrons, here 
pions, created in partial waves with a relative angular momentum difference of one. The 
dominant contribution being from the s-p interference term [10,11]. Therefore information 
about the IFF can be gained from a partial wave analysis for di-pion production [12]. Here 
the data from [13] suggest a phase shift around the p mass, leading to a sign change in the 
fragmentation function. In some models [14] the location of the phase shift around the 
mass of the p meson is caused by the interference of pion pairs produced in a p wave coming 
from the decay of the spin one p which interferes with the non-resonant background. Based 
on this model and estimation of particle yields from simulations, Bacchetta, Ceccopieri, 
Mukherjee and Radici [15] made model predictions for the magnitude of IFF Asymmetries 
at Belle. Again, a strong dependence on the invariant mass is predicted, with a maximum 
around the rho mass. The asymmetries are expected to rise with z due to the preservation 
of spin information early in the fragmentation. 



2.2 The Belle experiment 

The Belle detector, located at the KEKB asymmetric energy e~^e~ collider is described in 
detail in [16]. For the purpose of this study it is important that a high number of events 
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is recorded, a good particle identification allows to identify pions up to high values of z 
and that the detector is hermetic to minimize acceptance effects. Both is well fulfilled 
by the Belle detector. It is located at the collision point of the 3.5 GeV and 8 GeV 
e~ beam and is almost symmetric in the center of mass system of the beams. It is a 
large solid angle magnetic spectrometer consisting of barrel and end caps parts. For a 
more homogeneous acceptance function only the barrel part was used for this analysis. 
It comprises a silicon vertex detector, a 50-layer drift chamber and an electromagnetic 
calorimeter (CsTI) in a magnetic field of 1.5T. Particle identification over a wide range of 
momenta is done using an array of aerogel Cherenkov counters, a barrel-like arrangement 
of time-of-flight scintillation counters and an instrumented iron flux return yoke outside 
of the coil to detect K\ mesons and muons. 

2.3 Data selection and Asymmetry extraction 

From the 661 fb~^ data sample roughly 589 fb~^ were taken on the T(4S') resonance and 
73 fh~^ taken in the continuum 60 MeV below. Because the thrust cut used to select 
events with a two jet topology containing light and charm quarks also rejects events in 
which B mesons were produced, on-resonance and continuum data can be combined. The 

thrust is defined as T = ^ and fi is direction of the thrust chosen such that T is 

Ei \Pi\ 

maximal. B meson events produced on the T resonance have a spherical shape, since 
their high mass does not allow for large kinetic energy. On the other side, light quark anti 
quark pairs have a more two jet like topology.. An applied thrust cut of T > 0.8 reduces 
the contamination with B events to an order of 2% [5]. The inversion of the thrust cut 
selects events that don't have a clear two jet topology and that are contaminated by T 
decays into B mesons. This leads to a decrease of the asymmetry. 

As described earlier, only the barrel region of the detector is used in the analysis. 
Therefore the thrust axis is required to lie in a region well contained within it. This 
translates into a cut on the z component of the thrust axis \wz\ < 0.75, which also allows 
for the particles of the jet around the axis to be reconstructed in the barrel. Further cuts 
on the event level are a reconstructed energy of at least 7 GeV to reject e"'"e~ t^t~ 
and to reliably reconstruct the thrust axis. The later is computed using all charged tracks 
and photons passing some minimum energy cuts. A mean deviation of 135 mrad with a 
RMS of 90 mrad of the thrust axis from the real quark-anti quark axis is computed from 
simulations. 

Only events were selected which satisfy a vertex cut of 2 cm in the radial and 4 cm in 
the beam direction. On the track level the fiducial cuts to reduce acceptance effects are 
a constraint to the barrel region of the detector using a cut on the polar angle 6 in the 
laboratory system of —0.6 < cos(6') < 0.9 which translates to an almost symmetric cut in 
the CMS. To make sure that the azimuthal range of tracks around the thrust axis is not 
biased, only tracks are chosen that have at least 80% of their energy along the thrust axis. 
This restricts tracks to be within a cone that is entirely contained in the acceptance and 
reduces false asymmetries considerably as shown later in sec. 12. 4[ No false asymmetries 
from this cut is expected and since most of the energy of the jet is contained within the jet 
no significant dilution of the asymmetries either. Only tracks above a minimal fractional 
momentum z > 0.1 are considered. They have to be positively identified as pions. These 
tracks are then sorted into two hemispheres according to the sign of their momentum 
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projection on the thrust axis. All possible pairs of vr+vr" in the same hemisphere are 
selected and the angle </)q, a G {1, 2} computed. To this end the vector = Pi — P2 
for each pion pair in hemisphere a is formed and the azimuthal angle around the thrust 
axis with respect to the event plane computed according to eq. [H Charge ordering of 
/ii, h2 in the computation of R is always the same so that the effect is not averaged 
out. A weighting of the hadron momentum vector with the inverse fractional energy z as 
suggested by [17] only led to differences in the thrust axis within numerical uncertainties. 
For the computation of the yields in a specific 0i + 02 bin all combinations of pairs in 
the two hemispheres are considered. Due to the cone cut described earlier, the number 
of hadrons with a false hemisphere assignment is negligible. From Monte Carlo studies a 
signal purity for di-pion pairs (4 particles) of better than 90% over the whole kinematic 
range is obtained. 

2.4 Systematics studies 

In order to determine the systematic error on the measurement, simulation and real 
data was used to determine the contribution of detector effects and competing physical 
processes to false asymmetries or a dilution of the measurement. The studies that lead 
to the biggest contribution to the systematic error are the study of false asymmetries in 
simulation and real data in which no asymmetry is expected due to a wrong assignment of 
the thrust axis or the hemispheres of the particle pairs. Any false asymmetries, together 
with their statistical errors were added to the systematic error. 

Since the IFF effect is not included in the Pythia event generator used, checking for the 
asymmetries in fully simulated data allows to estimate effects of the detector acceptance 
and efficiencies on the asymmetry. Table [1] shows the false asymmetries extracted from a 
full simulation of the detector using GEANT. 

Another source for false asymmetries are events reconstructed from real data in which 
the asymmetries are averaged out. This is the case for mixed events in which the angles 0i 
and 02 are taken from different events or events in which the angels are computed for pion 
pairs in the same hemisphere. The later case leads to asymmetries due to phase space 
restrictions, which could be reproduced, the former to asymmetries in the order of one 
per-mille, shown in table [21 which has been added to the systematic error. Contributions 
from higher harmonics in the fits to the cosine modulation are under one per-mille. 

Due to the smearing of the thrust axis reconstruction with respect to the true quark- 
anti quark axis, the extracted asymmetries are diluted. Since this dilution could be 
reproduced in weighted Monte Carlo using the observed smearing of the thrust axis, it 
was corrected for. 

More physics motivated systematic checks concern the dependence on the various 
kinematic factors. The dependence on the kinematic factor ^^^'"^ J ^ should be linear, as 
should be the dependence on sin Q^^ if the effect is dominated by the s-p interference term. 
Both could be validated. 

Even though the Belle detector is very stable over time, checks have been done to 
determine the compatibility of data taking periods and data taken on and off the T 
resonance. To this end, the of each fit has been computed and all values have been fitted 
by an appropriate distribution. The result of these tests show very good compatibility. 

We checked for correlations in the data that might lead to false error estimates by 
breaking up weighted Monte Carlo data in 500 small chunks and comparing the mean 
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error of the extracted asymmetries to their variation. The results are compatible, so no 
systematic error was assigned. 

For the interpretation of the results the respective fraction of processes contributing 
to the asymmetries are very important. These have been estimated from Monte Carlo 
simulations and are shown in fig. [21 It is evident that the contribution from T decays has 
been almost eliminated by the thrust cut. There are also marginal contributions from 
r pairs. For these the false asymmetries are compatible with zero and they have been 
added with their statistical error to the overall systematic error. Besides the contribution 
from light quark pairs there is a considerable contribution from charm quarks. This 
contributions decreases with increasing z which can be understood, since charmed mesons 
have to undergo an additional decay before they can contribute to the pion asymmetries. 
This decay can also cause the invariant mass dependence namely a general decrease with 
higher masses. However, the highest bin shows a very high charm contribution, in some 
bins more than half of the events, which is not yet understood. The asymmetry results 
for these bins indicate that the IFF for charm quarks is non-vanishing and of similar 
magnitude as that of lighter quarks. 












(a) 



(b) 



Figure 2. Relative process contributions from light quark-anti quark events (red), charm events (green), 
charged B meson pairs (blue) , neutral B meson pairs (yellow) and r pairs (purple) as a function of zi for 
all z\ bins (a) and as a function of 7712 for all m \ bins (b) 



2.5 Results 

The results obtained for the a\2 asymmetry as defined in eq. [2]are shown in figs. [3] binned 
in the invariant masses mi, m2 of the hadrons pairs in the first and second hemisphere and 
their fractional energies zi , z^^ respectively. Table [3] shows the integrated asymmetries 
and the averaged kinematic observables. The extracted asymmetries are large, especially 
when considering that a product of the IFF for quark and anti quarks is measured. As 
expected the magnitude of the effect rises with z. However, the invariant mass behavior 
does not match model predictions from [15]. But these model predictions are only available 
in leading order and heavily dependent on simulations which were tuned for the SIDIS 
experiment HERMES at a center of mass energy of roughly 7 GeV. The asymmetries 
rise up to around the mass of the p but then plateau instead of decreasing again. A 
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sign change of the IFF can therefore not be confirmed. However bins with high invariant 
masses receive also considerable contributions from charm quarks as shown before. 
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Figure 3. Results for the ai2 modulations in a symmetric 8x8 binning in rrii, 7712 for TO{i_2} between 
0.25 GeV and 2 GcV (a) and for a symmetric 9x9 binning in zi, Z2 for Z{i,2} between 0.2 and 1 (b). 
The statistical error is shown in blue, the systematic error in green 
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2.6 Summary and Outlook 



The first direct measurement of tlie interference fragmentation function using 661 fb~^ of 
data recorded at the Belle experiment has been presented. The asymmetries are large, up 
to 10 %, which would correspond to an IFF contribution of over 30%. Models predicting 
a sign change or a decrease for invariant masses higher than the p meson's could not be 
confirmed. Charm quarks play a significant role at high invariant masses They seem to 
introduce an IFF asymmetry of similar magnitude as light quarks and further studies 
are under way to determine the charm quark contribution to the asymmetries. The 
analysis presented here should enable a combined analysis to extract the transversity 
distribution of data taken in SIDIS, proton-proton and e~^e~. This is very desirable 
due to the various advantages as compared with the extraction via the Collins effect 
and its complementarity. In proton-proton collisions the use of the IFF effect to access 
transversity is especially helpful, since it can help disentangle different contributions to 
the measured transverse single spin asymmetries Ajy. Our plans for the future contain 
also an extraction of other particle combinations, namely those including neutral pions 
and charged kaons. Furthermore an extraction of the unpolarized yields is planned to 
facilitate the extraction of the IFF. 



Table 1. MC results averaged over all z bins in %. 



sample 


species 


Zi , 2;2-Asymmetries 






(012) 


(ai2i?,) 




No 


opening cut 




uds 47r 


TTTT 


-0.089 ±0.008 


-0.108 ±0.008 


uds acceptance 


TTTT 


-0.488 ±0.011 


-0.490 ±0.011 


uds MC rec. 


TTTT 


-0.394 ±0.013 


-0.418 ±0.013 


charm rec. 


TTTT 


-0.446 ± 0.041 


-0.388 ± 0.044 


With opening cut of 0.8 


uds 4:71 


TTTT 


-0.038 ±0.013 


-0.035 ±0.013 


uds acceptance 


TTTT 


-0.112 ±0.016 


-0.113 ±0.016 


uds MC rec. 


TTTT 


0.012 ±0.019 


0.008 ±0.019 


charm rec. 


TTTT 


0.006 ± 0.040 


0.027 ±0.040 



Table 2. Mixing results averaged over the z binning in %. The results integrated over other binnings 
are nearly identical. 



sample 


zi , Z2-Asymmetries 








uds An 


0.070 ±0.013 


0.030 ±0.013 


uds acceptance 


0.020 ±0.016 


-0.021 ±0.016 


uds rec. 


0.091 ±0.019 


0.087 ±0.019 


charm rec. 


-0.024 ±0.040 


-0.017 ±0.040 


Data 


-0.019 ±0.017 


-0.012 ±0.017 
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Table 3. Integrated asymmetries and average kinematics. 



(^i), {Z2) 0.4313 

(mi), (ma) 0.6186 

(sin2^/(l + cos2e)) 0.7636 

(sin^;,) 0.9246 

(cos^,,) 0.0013 



-0.0199 ±0.0002 ±0.0009 



References 

J. C. Collins, S. F. Heppelmann and G. A. Ladinsky, Nucl. Phys. B 420, 565 (1994) 
[arXiv:hi^^/9305309j . 

A. Airapetian et al. [HERMES Collaboration], JHEP 0806 (2008) 017 
j arXiv:0803.2 367l [hep-ex]]. 

Talk given at DIS2009 by H. Wollny for the COMPASS collaboration 

Talk given at PKU-RBRC Workshop 2008 by R. Yang for the PHENIX collaboration 

R. Seidl et al. [Belle Collaboration], Phys. Rev. D 78, 032011 (2008) |arXiv:0805.2975] 
[hep-ex]]. 

M. Anselmino, M. Boglione, U. D'Alesio, A. Kotzinian, F. Murgia, A. Prokudin and 
S. Melis, arXiv:0812.4366l [hep-ph]. 

F. A. Ceccopieri, M. Radici and A. Bacchetta, Phys. Lett. B 650, 81 (2007) 
[arXiv:hep-ph/07032 65| . 
D. Boer, Nucl Phys B603, 192 (2001) 

D. Boer, R. Jakob and M. Radici, Phys. Rev. D 67 (2003) 094003 
[arXiv:hep-ph/0302232] . 

Radici, M. and Jakob,R. and Bianconi,A.,Phys. Rev.D65,074031(2002) 
Bianconi,A. et al. Phys. Rev.D62,034009(2000) 
Jaffe, R.L. and Jin, X. and Tang,J.,Phys. Rev. Lett. 80,1166, (1998) 
Estabrooks,P. and Martin,A.D.,Nucl. Phys.,B79,301(1974) 
Bacchetta,A. and Radici,M.,Phys. Rev.,D74,114007(2006) 

A. Bacchetta, F. A. Ceccopieri, A. Mukherjee and M. Radici, Phys. Rev. D 79, 
034029 (2009) | arXiv:0812.0 6TT] [hep-ph]]. 

A. Abashian et al. (Belle Collab.), Nucl. Instr. and Meth. A 479, 117 (2002). 
X. Artru and J. C. Collins, Z. Phys. C 69 (1996) 277 |arXiv:hep-ph/9504220| . 



9 



